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The possibility of multi-band conductivity and multi-gap superconductivity is explored in oriented 
' VsSi thin films by means of reflectance and transmittance measurements at terahertz frequencies. 

The temperature dependence of the transmittance spectra in the normal state gives evidence of two 
bands contributing to the film conductivity. This outcome is consistent with electronic structure 
calculations performed within density functional theory. On this basis, we performed a detailed 
data analysis and found that all optical data can be consistently accounted for within a two-band 
^ ■ framework, with the presence of two optical gaps in the superconducting state corresponding to 

5^ ' 2A/fcT c values close to 1.8 and 3.8. 

. PACS numbers: 74.70.Ad, 74.25.Gz, 78.30.-j 

Large interest has been devoted p], Q to multi-band superconductivity after the discovery of two bands with 
two distinct superconducting gaps in MgPj2 3]. Very recently, the interest in multi-band superconductivity has 
been renewed 0] by the discovery of the Fe-As based superconductors, where the presence of multiple bands is well 
established @. It is worth noting that two-gap superconductivity has been considered theoretically since the fifties 
@, and experimentally observed in transition metals 0, in the A15 compound NbaSn [1, Q and in MgCNi 3 Qjl . 
However, only after the discovery of superconductivity in MgB2, clear cut evidence of two-gap superconductivity and 
of its implications was obtained, since in this compound the multi-band, multi-gap character is emphasized by the 
! \ exceptionally low inter-band scattering [ill ]. 

In the case of the A15 VsSi system, a spread in the gap values corresponding to 2A/fcT c extending from 1.0 to 3.9 
' has been reported in the past [H, [l3| . Recently, the electrodynamic response in the microwave region gave evidence 
C*~) . of two gaps [HI , while muon spin rotation measurements were consistent with a single-gap model [lj| . It is also worth 
• noting that VsSi has been recently treated as a two-band system [l6l [l7|. but existing band structure calculations 
[let fig ] do not provide detailed enough information. 
[ — ■ In this contradictory scenario, we investigate in the present work the possible two-band, two-gap character of V^Si, 
by performing both an infrared spectroscopy study and electronic structure calculations within density functional 
theory (DFT). 

Infrared (IR) spectroscopy is a powerful tool to study the properties of a conducting system. In the normal state 
. (N-state), the Drude model for the frequency-dependent complex conductivity ctjv = <J\n + i&2N can be employed to 
describe the optical response of free-charge carriers [2(| • As shown in the MgPj2 case [U [22| , it is thus possible to 
determine the contributions of different i-bands and thus the corresponding plasma frequencies tti and scattering rates 
7i. In the superconducting state (S-state), on the basis of the BCS model for the complex conductivity as = <Jis+i<J2S 
(231 ] . far-IR/terahertz measurements can be of particular importance since a mark of the superconducting gap A can 
be observed at Tilu ~ 2A (optical gap) for an isotropic s-wave BCS superconductor. In particular, a maximum at the 
optical gap is expected either in the ratio Rs/Rn (for a bulk sample) or in the ratio Ts/Tn (for a thin film), where 
Rs (7s) and Rn (Xn) are the frequency-dependent reflectances (transmittances) in the S- and N- state, respectively 
[23 ]. In the MgB 2 case [2l|, evidence of the two gaps was detected in the Rs/Rn spectrum of an ultra-clean film [22| 
by assuming a parallel sum of the conductivity of two independent bands and by using the model cr<j introduced by 
Zimmermann et al. [24j which generalizes the BCS one to arbitrary T and 7 values (Zimmermann model). 

We performed IR measurements on high-quality V^Si textured films. Details on the film growth by pulsed laser 
deposition and on their properties are reported elsewhere [25|. We studied two films grown on LaAlOs (LAO) (001) 
0.5 mm thick substrates, which exhibit preferential (210) orientation along the out-of-plane direction. The first film 
of thickness d=180 nm (film dl80) has good transport properties (resistivity at 300 K close to 200 /ifkm, residual 
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FIG. 1: (Color online) (a) Transmittance spectra T(w) of film d033 at selected temperatures in the THz region. The 20 K 
spectrum is compared with the two-band (2-b) best fit. (b) T(uj) spectra of a LaA103 substrate (0.5 mm thick) at the same 
temperatures [261 ]. 



resistivity ratio RRR=8) and T c =16.1 K; the second film, 33 nm thick (film d033), has worst transport properties 
{RRR=A.h) and a slightly lower T c value (T c =15.3 K). 

We first performed measurements of the IR reflectance R(w) of film dl80 as a function of temperature (not shown). 
The R(ui) spectrum at 300 K is in very good agreement with the result of previous measurements We attempted 
an analysis of the measured spectra by including in the model ffjv a number of Lorentz contributions representing 
optically active inter-band transition. However, unambiguous results were not obtained since at any temperature a 
very broad Drude contribution [19j is partially overwhelmed by strong interband transitions occurring at rather low 
frequencies because of the complex electronic structure of VaSi [l8L Ha]. 



We then measured the transmittance T(u>) spectrum of film d033 by focusing the IR beam on the film or on a 
hole for reference. Measurements were performed between 300 and 20 K in the terahertz (THz) region (here defined 
by photon energies below 16 meV, frequency to < 130 cm -1 or v < 4 THz), where the LAO substrate is partially 
transparent. T{oS) spectra of film d033 are reported in Fig. QJi at selected temperatures, those of LAO [26| at the same 
temperatures in FigQJ) for comparison. The spectral shape of the 300 K T{u>) spectrum is qualitatively similar to that 
of LAO as it monotonically increases with decreasing frequency, while the transmitted intensity is strongly reduced 
by the above discussed broad Drude contribution in the film conductivity. On the contrary, the spectral shape of the 
20 K T(u>) shows a broad maximum and then decreases with decreasing frequency. This result can be qualitatively 
explained by the presence of a second sharp Drude contribution, which narrows with decreasing temperature thus 
reducing the transmitted intensity at very low frequencies. 

In order to better understand this result, we carried out electronic structure calculations within DFT on V3SL We 
calculated the Fermi velocities starting from the energy bands shown in the inset of FigO We used the experimental 
lattice constants a and c [27|, and calculated the dimerization of V chains along x and y directions to be 0.0025a. 
The Fermi level Ep of VsSi is crossed by four bands (two bands with similar character and small density of states are 
grouped together) with mostly V 3d character; with the chain direction along z, orbitals d x 2_ y 2 and d z i on one side 
and d xz and d yz on the other produce narrow and wide bands, respectively. 

In Fig. [2] we plot the distribution of values of the moduli of Fermi velocities [28[ for the bands crossing Ep; the 
distributions are normalized to the densities of states of the same bands, so that the average squared velocity times the 
integral of the distribution provides, apart from numerical factors, the contribution to the square of the total plasma 
frequency Sltot = 3.34 eV. While it is difficult to assign a definite orbital character to each band throughout the whole 
Brillouin zone, we can see that the band separation naturally provides different average velocities, showing an overall 
bimodal distribution given by the bands n. 3 and 4, while the remaining bands provide a smaller contribution to 
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FIG. 2: (Color online) Distribution of moduli of the Fermi velocities for the bands crossing Ep. Inset: Energy bands along the 
symmetry lines of the Brillouin zone; the different symbols identify the bands as in the inset. 



plasma frequencies due to their smaller density of states. The bimodal distribution resulting from our calculations 
indicates that the V^Si conductivity in the N-state can be safely described in terms of two bands characterized by 
different f2 values (two-band, 2-b model). It remains to be verified whether the presence of the two bands reflects in 
the opening of two distinct superconducting gaps, as in the MgB2 case. 



To this aim we measured at T < 20 K the R S (T)/R N (where R N = R(20 K)) spectra of the dI80 film (see Figdi,) 
and the Ts(T)/Tn (where T/v = T(20 K)) spectra of the d033 film (see FigHK). These measurements were made by 
cycling the temperature in the 6-20 K range, without collecting reference spectra. In this way one avoids any variation 
in the sample position and orientation, which may yield frequency-dependent systematic errors in R(uj) and T(uj) [29| . 
All measurements in the THz region were made by employing synchrotron radiation at the infrared beamline SISSI 
[30l | at the synchrotron Elettra (Trieste, Italy) and at the BESSY storage ring (Berlin, Germany), where coherent 
synchrotron radiation is available [3l| . We remark that a high-flux synchrotron source allows an high accuracy in the 
detection of small effects in the R(T)/Rn spectra, and overcomes the problem of the very low intensity transmitted by 
the film+substrate system in the case of transmission measurements on a conducting film. We also note that 7s /7/v 
(where 7s = 7~(6 K)) is affected by the superconducting transition much more than Rs/Rn (where Rs — R(Q K)) 
since the transmitted intensity is dominated by the effect of the absorptive processes in the film, which significantly 
changes below T c [H,^. 



We first notice that the R{T)/Rn spectrum increases on decreasing temperature (see Fig. [3h), until Rs/Rn reaches 
a maximum and becomes nearly constant below 20 cm . This indicates the presence of a superconducting gap A 
close to 10 cm -1 . Indeed, for u> — > 0, the reflectance Rn of a conducting system tends to 1 for a bulk system, to a 
slightly lower value for a thin film. Therefore, since Rs approaches 1 at u> — 2 A, Rs/Rn exhibits a maximum around 
2A in the case of a bulk sample, remains nearly constant below 2A in the film case. As to the T(T)/Tn data (see 
Fig. H^,), a maximum develops on decreasing T until the 7s /7/v exhibits a well defined peak around 40 cm -1 , which 
indicates the presence of a superconducting gap A around 20 cm . 

For a detailed analysis of the measured spectra based on the 2-b model, the conductivity of VaSi can be described 
by the sum of two Drude terms in the N-state (oat, with parameters OjjT'j, i = A,B), and by the sum of two 
Zimmermann terms in the S-state, as discussed above (as, with parameters T, f2j, 7^, Aj, i = A, B). From cr, by using 
standard relations (2(|, it is possible to compute the model refractive index n = n + ik of V^Si in both the N- and 
S-state. The transmittance and reflectance spectra of the film+substrate system in both states can then be evaluated 
by means of an exact procedure [33| which requires, besides thickness, numerical values of n and k for both film and 
substrate. For LAO, n and k values in the THz region were obtained from previous measurements of transmittance 
and reflectance of a LAO substrate [26| through a numerical, model independent procedure [34| . 
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FIG. 3: (Color online) (a) R(T)/Rn spectra of film dl80 at selected temperatures in the Thz region, (b) Rs/Rn (with 
Rs = R(6K)) spectrum compared with the two-band (2-b) and one-band (1-b) best fits. In the I-b case, best-fit parameters 
are SI = 2.35 eV, 7 = 125 cm -1 , and A = 14 cm -1 , (c) ais (in units 10 4 f2 _1 cm _1 ) of VaSi from 2-b model. 
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FIG. 4: (Color online) (a) T(T)/T N spectra (with T N =T(2Q K)) of film d033 at selected temperatures in the THz region, (b) 
Ts/Tn spectrum compared with the two-band (2-b) best fit (SI a = 3.35 eV, 7a= 1000 cm -1 , Aa~ 12 cm -1 , SIb = 1-24 eV, 
jb~ 75 cm -1 , Ab~ 21 cm -1 ), (c) absorption coefficient as (in units 10 4 cm -1 ) of VaSi from 2-b model. 
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We first analyzed the Rs/Rn spectrum, by performing a 6 parameter fit (f^,7i, Aj, i = A, B). The best-fit curve 
does well describe experimental data, as shown in FigJSJa. By keeping into account the uncertainties resulting from 
the fitting procedure, we can safely pose £l A = 3.7 ± 0.1 eV, D, B = 1.00 ± 0.05 eV, j A = 800 ± 40 cm" 1 , lB = 45 ± 4 
cm -1 , = 10±1 cm -1 , and Ab — 19±2cm _1 . We verified that unsatisfactory results are obtained by considering 
one band only (one-band, 1-b model). Note that, in the 2-b approach, the A-band contribution has a crucial effect, 
as it gives a vanishing a±s only below 2Aa (see Figj3fc). Since, in general, Rs approaches 1 when a\s vanishes, this 
explains why only the gap is well evident in the Rs/Rn spectrum. 

In order to verify the compatibility of this fit with first principle results, we calculated plasma frequencies from our 
Fermi velocities. In particular, we integrated over the Brillouin zone separating the low and high vp regions (setting 
an arbitrary but plausible cutoff at 10 7 cm/s; see Fig(2]). While this procedure is only qualitative, it may be justified 
by the presence of clearly separable ranges of Vp, in turn deriving from different orbital natures. We obtain £1^=3.21 
eV and flB = 0.81 cV. In comparing these results with those of the fitting procedure, it is worth noting that the best-fit 
values of Ha and Ob, and thus of Sltot, can be overestimated because the interband transitions discussed above may 
give non negligible contributions to &n even at very low frequencies. This effect is minimized in the best-fit Qa/Qb 
ratio (3.7±0.3), which results to be in a remarkably good agreement with the computed one (3.96). 

In analysing the transmittance spectra of film d033, we verified that only the 2-b model well describes T{oj) at all 
temperatures, but the resulting £li and ji values are not unambiguously determined. We thus used a procedure which 
simultaneously fits both Tn and Ts/Tn, thus imposing important constraints in the fitting procedure. Good fits of 
both TJv (see Fig0J>) and Ts/Tn (see FigHja) were obtained, with f2j values slightly lower, 7, values slightly higher 
than those found in fitting the Rs/Rn spectrum. This can be simply explained by the worse conducting properties 
of film d033 with respect to film dl80. As to the gap values, nearly equivalent fits are obtained for A^ ranging 
from 11 to 16 cm" 1 , while Ab = 21.0 ±0.5 cm -1 is well determined since it corresponds to the peak in Ts/Tn- The 
transmission measurement, dominated by the absorption process, thus permits to unambiguously establish the Ab 
value which was more poorly defined in the Rs/Rn measurement. The higher sensitivity of Ts/Tn to the Ab gap 
is a consequence of the frequency dep endence of the absorption coefficient as of VsSi in the S-state. Indeed, as, as 
evaluated with standard relations [201 ] through the 2-b model, only exhibits a well defined minimum around 2A^ (see 
Figgb). 

In summary, we addressed the debated problem of multi-band, multi-gap nature of VaSi by means of reflectance 
and transmittance measurements in the THz region on high quality oriented films. Experimental results indicate 
the presence of two bands contributing to the VaSi conductivity in the normal state, and of two optical gaps in the 
superconducting state. Electronic structure calculations within density functional theory showed that the distribution 
of the modulus of the Fermi velocity exhibits a clear bimodal character, which indicates that the VsSi conductivity 
in the normal state can be safely described in terms of two bands, characterized by different plasma frequencies. On 
this basis, we performed a detailed data analysis and found that all optical data can be consistently accounted for 
within a two-band framework, with the presence of two optical gaps corresponding to 2A/kT c values close to 1.8 and 
3.8. 

The authors acknowledge M. Prasciolu for preparing the Au reference surface used in the absolute reflectance 
measurements. 



[1] F. Bouquet, Y. Wang, I. Sheikin, T. Plackowski, A. Junod, S. Lee, and S. Tajima, Phys. Rev. Lett. 89, 257001 (2002) 
[2] I.I. Mazin and VP. Antropov, Physica C 385, 49 (2003) 

[3] J. Kortus, I.I. Mazin, K.D. Belashchenko, VP. Antrop ov, and L.L. Boye r, Phys. Rev. Lett. 86, 4656 (2001) 

[4] A. Moreo, M. Daghofer, A. Nicholson, and E. Dagotto. larXiv:0906. 17471 

[5] I.I. Mazin, D.J. Singh, M.D. Johannes, and M.H. Du, Phys. Rev. Lett. 101, 057003 (2008) 

[6] H. Suhl, B.T. Matthias, and L.R. Walker, Phys. Rev. Lett. 3, 552 (1959). 

[7] L.J. Vieland and A.W. Wicklund, Phys. Rev. 166, 424 (1968). 

[8] J.C.F. Brock, Solid State Commun. 7, 1789 (1969) 

[9] V. Guritanu, W. Goldacker, F. Bouquet, Y. Wang, R. Lofiz, G. Goll and A. Junod, Phys. Rev. B 70, 184526 (2004) 
[10] A. Walte, G. Fuchs, K.-H. Muller, A. Handstein, K. Nenkov, V.N. Narozhnyi, S.-L. Drechsler, S. Shulga, L. Schultz, and 

H. Rosner, Phys. Rev. B 70, 174503 (2004). 
[11] I.I. Mazin, O.K. Andersen, O. Jepsen, O.V. Dolgov, J. Kortus, A. A. Golubov, A.B. Kuzmenko, and D. Van Der Marel, 

Phys. Rev. Lett. 89, 107002, (2002). 
[12] D.B. Tanner and A.J. Sievers, Phys. Rev. B 8, 1978 (1973) 

[13] For a number of gap values measured in VaSi see: A.B. Mitrovic and JP. Carbotte, Phys. Rev. B 26, 1244 (1982) 

[14] YA. Nefyodov, A.M. Shuvaev, and M.R. Trunin, Europhys. Lett. 72, 638 (2005). 

[15] F.D. Callaghan, M. Laulajainen, CV. Kaiser, and J.E. Sonier, Phys. Rev. Lett. 95, 197001 (2005). 

[16] N. Hauptmann, M. Becker, J. Kroger, R. Berndt, Phys. Rev. B 79, 144522 (2009) 



6 



[17] V.G.Kogan, C.Martin, and R.Prozorov, arXiv:0905.0029 

[18] B.M. Klein, L.L. Boyer, D.A. Papaconstantopoulos, and L.F. Mattheiss, Phys. Rev. B 18, 6411 (1978). 
[19] A. Borghesi, A. Piaggi, G. Guazzetti, F. Nava, and M. Bacchetti, Phys. Rev. B 40, 3249 (1989). 
[20] G. Burns, Solid state physics, Academic Press, Boston (1990). 

[21] For a comprehensive review of infrared studies on MgB2, see: A.B. Kuzmenko, Phys. C 456, 63 (2007). 

[22] M. Ortolani, P. Dore, D. Di Castro, A. Perucchi, S. Lupi, V.Ferrando, M.Putti, I. Pallecchi, C. Ferdeghini, and X.X.Xi, 

Phys. Rev. B 77, 100507(R) (2008). 
[23] M. Tinkham, Introduction to Superconductivity, McGraw-Hill, New York (1975). 
[24] W. Zimmermann, E.H. Brandt, M. Bauer, E. Seider, and L. Genzel, Phys. C 183, 99 (1991). 

[25] C. Ferdeghini, E. Bellingeri, C. Fanciulli, M. Ferretti, P. Manfrinetti, I. Pallecchi, M. Putti, C. Tarantini, M. Tropeano, A. 

Andreone, G. Lamura, and R. Vaglio, IEEE Trans. Appl. Supercond. (in press, 2009). 
[26] P. Dore, G.P. Gallerano, A. Doria, E. Giovenale, R. Trippetti, and V. Boffa, Nuovo Cimento D 16, 1803 (1994). 
[27] P. Chaddah and R.O. Simmons, Phys. Rev. B 27, 119 (1983) 

[28] A similar distribution in terms of the components of vf would peak around zero due to geometrical factors. 

[29] M. Ortolani, S. Lupi, L. Baldassarre, U. Schade, P. Calvani, Y. Takano, M. Nagao, T. Takenouchi, and H. Kawarada, 

Phys. Rev. Lett. 97, 097002 (2006) 
[30] S. Lupi, A. Nucara, A. Perucchi, P. Calvani, M. Ortolani, L. Quaroni, and M. Kiskinova.JOSA B, 24, 959 (2007) 
[31] W.B. Peatman and U. Schade, Rev. Sci. Instr. 72, 1620 (2001) 

[32] G.P. Williams, R.C. Budhani, C.J. Hirschmugl, G.L. Carr, S. Perkowitz, B. Lou, and T.R. Young, Phys. Rev. B 41, 4752 
(1990) 

[33] P. Berberich, M. Chiusuri, S. Cunsolo, P. Dore, H. Kinder, and CP. Varsamis, Infrared Phys. 34, 269 (1993). 
[34] S. Cunsolo, P. Dore, and CP. Varsamis, Applied Optics 31, 4554 (1992) 



